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Abstract 



The ratio of neutral Higgs field vacuum expectation values, tan/3, is one of the most important 
parameters to determine in either the Minimal Supersymmetric Standard Model (MSSM) or a general 
type-II Two-Higgs Doublet Model (2HDM). Assuming an energy and integrated luminosity of y/s = 
500 GeV and C = 2000 fb"^ at a future hnear colUder (LC), we show that a very accurate determination 
of tan P will be possible for low and high tan (3 values by measuring the production rates of Higgs 
bosons and reconstructing Higgs boson decays. In particular, based on a TESLA simulation, and 
assuming no other light Higgs bosons and 100 < niA < 200 GeV, we find that the rate for the process 
e^e~ — > bbA — > bbbb provides a good determination of tan/3 at high tan f3. In the MSSM Higgs sector, 
in the sample case of = 200 GeV, we find that the rates for e~^e~ bbA + bbH bbbb and for 
e^e~ — > HA bbbb provide a good determination of tan /3 at high and low tan /?, respectively. We also 
show that the direct measurement of the average total widths of the H and A in e~^e~ — > HA — > bbbb 
events provides an excellent determination of tan/3 at large values. In addition, the charged Higgs 
boson process e~^e~ H~^H~ — > tbtb has been studied. The sensitivity to tan /3 at the LHC obtained 
directly from heavy Higgs boson production is briefly compared to the LC results. 
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I. INTRODUCTION 

Theories beyond the Standard Model (SM) that resolve the hierarchy and fine-tuning prob- 
lems typically involve extensions of its single- doublet Higgs sector to at least a two-doublet 
Higgs sector (2HDM) The most attractive such model is the Minimal Supersymmetric 
Standard Model (MSSM), which contains a constrained two-Higgs-doublet sector 0. In other 
cases, the effective theory below some energy scale is equivalent to a 2HDM extension of the 
SM with no other new physics. Searching for the Higgs particles and studying their properties 
have high priority for both theoretical and experimental activities in high energy physics. 

Among other new parameters in 2HDM and SUSY theories, one is of particular importance: 
the ratio of the vacuum expectation values of the two Higgs fields, commonly denoted as 
tan/3 = V2/V1. It characterizes the relative fraction that the two Higgs doublets contribute to 
the electroweak symmetry breaking = v1 + v^-, where v ~ 246 GeV. The five physical Higgs 
states couple to the fermions at tree-level as 

.mtcosa mt j .mbsina .mb 
htt : —I ^ —I — hho : % ^ —% — (1) 

f sin p t> V cos p V 

.m^sina .rrit ^ .mucosa .rrih ^ 
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Alt : -— cot 13 75 Ahh : -— tan 13 75 (3) 

V V 

H+ib : i-^Ktan/5(l + 75) +mtCot/3(l -75)], (4) 
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where a is the mixing angle in the CP-even sector, and the approximation indicates the decou- 
pling limit for itla ^ ttlz in the MSSM p, |], in which the couplings of the light Higgs boson 
h become SM-like. Eqs. show that tan/3 governs the coupling strength of Yukawa 

interactions between the fermions and the heavy Higgs bosons. In fact, heavy Higgs boson 
measurements sensitive to their Yukawa couplings are far and away the most direct way to 
probe the structure of the vacuum state of the model as characterized by the ratio of vacuum 
expectation values that defines tan /3. 

The parameter tan/3 enters all other sectors of the theory in a less direct way [|I|. For 
instance, in supersymmetric theories the interactions of the SUSY particles have tan j3 depen- 
dence. In addition, the relations of SUSY particle masses to the soft SUSY breaking parameters 
of supersymmetry involve tan (3. The renormalization group evolution of the Yukawa couplings 
from the unification scale to the electroweak scale is sensitive to the value of tan /3. The large 
top quark mass can be naturally explained with — rrir unification as a quasi-infrared fixed 
point of the top Yukawa coupling if tan/3 ~ 2 or tan/3 ~ 56 [^. The possibility of SO(IO) 
Yukawa unification requires high tan /? solutions . The predicted mass of the lightest SUSY 
Higgs boson also depends on tan (3, with a higher mass at larger tan /3 0. It will be very impor- 
tant to compare the measurements of and constraints on tan /3 from these other sectors of the 
theory to the direct determination of tan (3 coming from the heavy Higgs boson measurements 
that depend fundamentally on tan /? through the Yukawa couplings. 

Currently, some regions of the MSSM parameter space have been excluded at LEP due to 
the lower bound on the lightest Higgs boson mass. (A review of LEP-1 Higgs results shows 
possible signatures for all neutral and charged Higgs boson search channels [H].) Particularly 
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interesting is the exclusion 0.5 < tan (3 < 2.4 in the maximal-top-squark-mixing scenario when 
the SUSY scale Msusy ^ 1 TeV 0. More general MSSM parameter scans reduce the excluded 
tan/? range [|10|, especially if the top quark mass is allowed to vary within its la error range. 
Searches for top decay t — >■ H~^b at the Fermilab Tevatron, sensitive to tan/3 > 50 with 
rrit > friHi, and searches for the final state bbh bbbb, sensitive to tan/5 > 35, have also set 
limits [ITT] on the very high values of tan /3 as a function of the Higgs mass. Precision electroweak 
measurements may provide some weak constraints on tan (3 . Much of the parameter space, 
however, remains to be explored at future collider experiments. 

Because of the significance of tan j3 for the theory and phenomenology, it is important 
to constrain it and eventually to determine its value in future collider experiments. At the 
upgraded Tevatron (Run II) with high luminosity (10 fb~^ or higher), complementarity among 
the processes qq — > Wh,Zh |jl3|, gg — >• bbh,bbH,bbA |l5l, and gg — > h,H,A — > r+r" 



may allow SUSY Higgs detection throughout the full SUSY parameter (m^,tan/3) plane. 
Depending upon the integrated luminosity and the value of tan/3, either we will be able to 
directly observe the heavy Higgs processes and be able to determine tan j3 or, if the H, A, 
are not detected, we will be able to place an upper bound on tan/3 as a function of m^. At 
high tan/3, SUSY particle production may provide an additional handle [0 for exploration 
of the parameters as well. More recently, it has been pointed out that a large value of tan /3 
can substantially enhance B meson decay branching fractions |]18| and thus could enhance our 
ability to probe SUSY parameters through indirect SUSY and Higgs signals [|19]. SUSY Higgs 



detection at the LHC has been studied for many years; see for a summary of the early work. 
Recent studies by the LHC CMS and ATLAS collaborations can be found in |21|, |2^, |23|. The 
conclusion of these studies is that there is a "no-lose" theorem for SUSY Higgs discovery at 
the LHC, although there is a substantial region of (m^,tan/3) parameter space where only the 
light CP-even h will be detected. The determination of tan (3 at the LHC by measuring Higgs 
boson production rates was analyzed in and has been explored in greater detail in p^] . 
The possibility of measuring tan /3 at the LHC via production of gauginos and squarks was 
explored in |^, A future linear collider has great potential for discovering new particles 
and measuring their properties due to its clean experimental environment. An early study 
27| has been made to determine tan/? via gaugino production in e~^e~ collisions, followed by 
a discussion for ej collisions In fact, the value of tan/3 can be analytically determined 

by measuring parameters in the gaugino sector of the MSSM, as outlined in ||29|. Due to the 
clean experimental environment, stau pair production in e~^e~ collisions can be exploited to 
probe the SUSY parameters, in particular tan/?, via left-right mixing 0]. Especially useful at 
high tan /3 will be measurements of the r polarization in f decays, which is directly sensitive 
to mT-/(cos/3mvy) 

Many SUSY particle production and decay processes depend on sin /3 or cos j3 (and not their 
inverses). Once tan/? > 5, sin/3 ~ 1 and cos/3 is simply small, and thus even large changes in 
tan/3 will have little impact on the related experimental observables at large tan/? [2S]. There 
are exceptions to this statement, including the r polarization measurement referenced above. 
However, even this process involves other SUSY parameters at tree-level simultaneously as well. 
Thus, while many of the SUSY parameters can be measured with high precision, the accuracy 
with which tan j3 can be determined using measurements not involving the Higgs bosons remains 
uncertain, especially if tan /? is large. Further, in non-SUSY extensions of the SM, it may be 
that the only direct probe of tan/? will be via the Higgs sector. As seen from Eqs. (|I|)-(^), the 
heavy Higgs boson couplings to fermions are very sensitive to, and provide a direct probe of, 
tan /? at tree-level. Pair production of heavy Higgs bosons AH and H'^H~ was studied in e~^e~ 
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collisions [0, |3^, |3^ , and improved sensitivity to tan (3 was obtained by considering H^tb ^5 



and A/Hbb, A/Htt ||36|. If a muon collider becomes available to produce a heavy Higgs boson 
in the s-channel, its coupling could be measured with very excellent precision |^ . 

In this paper, we perform a comprehensive analysis of tan (3 determination via heavy Higgs 
boson production and decay at an e~^e~ linear collider with -y/i = 500 GeV and an integrated 
luminosity of 2000 fb~^. We amplify upon the results for the heavy neutral Higgs bosons 
obtained during the last Snowmass workshop (as summarized in |Q) and extend our study to 



include the charged Higgs boson. We show how various Higgs boson measurements can be used 
to determine tan (3. The different types of measurements we consider are complementary in that 
some provide good precision at low tan/3 and others at high tan/?; combined, a good determi- 
nation of tan j3 is possible throughout its whole range. We include background simulations and 
realistic 6-tagging efficiencies. In Sec. II, we focus on the heavy Higgs bremsstrahlung process 
bbA — > bbbb, the production rate for which is directly proportional to tan^ f3. We then include 
the bbH bbbb process in the context of the MSSM and estimate the accuracy with which 
tan (3 can be determined by combining experimental results for both processes. In Sec. HI, we 
examine the pair production of a CP-even Higgs boson [h or H) and the CP-odd A, followed by 
decay of the Higgs bosons to the bbbb final state. In particular, at large tan j3, the total decay 
widths of the heavy Higgs bosons can be broad since these widths are proportional to tan^ f3. 
The resulting accuracy for the tan/5 determination is obtained. We extend these studies to 
charged Higgs in Sec. IV. In Sec. V, we briefiy summarize the LHC sensitivity to tan (3 deriving 
from heavy Higgs production. Finally, we summarize our results in Sec. VI. 

Before proceeding with our analysis, we would like to point out that we are taking a phe- 
nomenological approach to the tan/? determination. Namely, we only consider tan/5 as an 
effective way of specifying the coupling for the Higgs bosons and fermions through the usual 
tree-level relations and explore the extent to which this coupling can be experimentally deter- 
mined at the linear collider experiments. After including radiative corrections, the relation of 
tan/? to the various Yukawa couplings becomes process-dependent. For a recent theoretical 
discussion on the issue of the gauge dependence of these relations, see Ref. In this context, 
our results should be viewed as giving the accuracy with which the actual Yukawa couplings 
can be measured. 

II. THE bbA bbbb BREMSSTRAHLUNG PROCESS 

Searches for bbA and bbh were performed in the four-jet channel using LEP data taken 



at the Z resonance In this section, we consider a linear collider with a 

center-of-mass energy of 500 GeV or higher, and begin by focusing just on the bbA pro- 
duction process that probes the direct coupling of the CP-odd Higgs A to bb. Our analysis 
will employ cuts designed to eliminate Higgs pair resonant production, which, when kine- 
matically accessible, dominates A production before cuts but is less sensitive to tan/?. The 
challenge of this study is the low expected production rate and the large irreducible back- 
ground for a four-jet final state, as discussed in a previous study |Q. A LC analysis has 



been performed using event generators for the signal process e~^e^ bbA bbbb and the 
e"'"e~ — >• eWu, e^e~Z, WW, ZZ, qq (g = u,d, s,c,b), tt, hA background processes |^ that 
include initial- state radiation and beamstrahlung. 

For a 100 GeV CP-odd Higgs boson and tan (3 = 50, the signal cross section is about 2 fb 
|8|, The generated events were passed through the fast detector simulation SGV [pO|]. 



The detector properties closely follow the TESLA detector Conceptual Design Report 
The simulation of the 6-tagging performance is very important for this analysis. The efficiency 
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FIG. 1: Simulated 6-tagging performance at the TESLA. 



versus purity distribution for the simulated b-tagging performance is shown in Fig. |l] for the 
hadronic event sample e'^e" — > qq for 5 flavors, where efficiency is the ratio of simulated bb 
events with the selection cuts to all simulated bb events, and purity is the ratio of simulated 
bb events with the selection cuts to all selected qq events. Details of the event selection and 



background reduction are described elsewhere [Q] 



For niA = 100 GeV in the context of the MSSM, the SM-like Higgs boson is the H while the 
light h is decoupled from WW, ZZ [cos(/9 — a) ~ 1 and sin(/5 — a) ~ 0]. The bbh coupling is 
essentially equal (in magnitude) to the bbA coupling (oc tan/5 at the tree level) and irih ~ itia, 
implying that the signal would be doubled from bbA and bbh. Also important will be hA pair 
production, which is proportional to cos(/3 — a) and will have full strength in this particular 
situation; HA production will be strongly suppressed. We focus first on bbA bbbb. 

The expected background rate for a given bbA bbbb signal efficiency is shown in Fig. ^ 
One component of the background is hA — >■ bbbb since it has rather weak dependence on tan/3. 
Our selection procedures are, in part, designed to reduce this piece of the background as much 
as possible. Nonetheless, it may lead to significant systematic error in the determination of 
tan/5 due to interference with the signal, as discussed below. For the bbA bbbb signal, the 
sensitivity S/ ^J~B for rrtA = 100 GeV is almost independent of the working point choice of signal 
efficiency in the range esei = 5% to 50%. For a working point choice of 10% efficiency, the total 
simulated background of about 16 million events is reduced to 100 background events with an 
equal number of signal events at tan (3 = 50. We estimate the error on determining tan (3 by 

A tan^ (3/ tan^ (3 = AS/S = ^fsTB/S. (5) 

If this were the only contributing process, then we would have \/ S + B / S ~ 0.14, resulting in 
an error on tan/3 = 50 of 7%. For smaller values of tan/3, the sensitivity decreases rapidly. A 
5a signal detection is still possible for tan/3 = 35. In the MSSM context, the bbh signal would 
essentially double the number of signal events and have exactly the same tan f3 dependence, 
yielding Atan^/^/tan^/? ~ v/300/200 ~ 0.085 for tan/3 = 50. 

Although the number of hA background events is very small compared to the other back- 
ground reactions after the event selection, interference between the signal {bbA — *■ bbbb plus 
bbh —>■ bbbb) and the background {hA — > bbbb) could be important. At the working point of 
10% signal efficiency, and after applying the selection procedures, the expected rate for the 
latter is 2 ± 1 events for C = 500 fb^^. To assess the effect of the interference, let us momen- 
tarily retain only the bbA signal and the hA background. We first calculate the cross sections 
cr(e+e" bbA — > bbbb), (T(e+e~ ^ hA ^ bbbb), and cr(e+e" — > bbA + hA ^ bbbb) with 
CompHEP 1 52] before selections. We define the interference as 

0"interf = (^bbA+hA ~ ^bbA ~ ^hA- (6) 
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For the default value = 4.62 GeV, at tan/? = 50 we obtain cr^^^ = 1.83 ± 0.01 fb, ahA = 
36.85 ± 0.10 fb, (Tf^iA+hA = 39.23 ± 0.12 fb, and thus dinterf = 0.55 ± 0.16 fb. We observe a 
constructive interference similar in size to the signal. Thus, more signal events are expected than 
simulated and the statistical error estimate is conservative. After selection cuts, we have found 
100 signal events versus two hA background events. The maximum interference magnitude 
arises if the interference events are signal-like, yielding an interference excess of (10 + ^/2Y - 
100 — 2 ~ 28, a percentage (~ 30%) similar to the ratio obtained before selection cuts. If the 
events from the interference are background-like, the resulting systematic error will be small, 
since the hA background is only a small part of the total background. In the MSSM context 
we have an exact prediction as a function of tan (3 for the combined contribution of hA bbbb 
and bbA — > bbbb (plus bbh — >■ bbbb) , including all interferences, and this exact prediction can be 
compared to the data. In order to test this exact prediction, it may be helpful to compare theory 
and experiment for several different event selection procedures, including ones that give more 
emphasis to the hA process. Of course, this exact prediction depends somewhat on other MSSM 
parameters, especially if decays of the h or H to pairs of supersymmetric particles are allowed 
or ratios of certain MSSM parameters are relatively large If this type of uncertainty 

exists, the systematic error on tan (3 can still be controlled by simultaneously simulating all 
sources of bbbb events for various tan/? values and fitting the complete data set (assuming 
that the other MSSM parameters are known sufficiently well). Another possible theoretical 
systematic uncertainty derives from higher-order corrections. The full next-to-leading order 



(NLO) QCD corrections are given in [gj, |5^. There it is found that using the running 6-quark 
mass incorporates the bulk of the NLO corrections. For example, for = 100 GeV, employing 
m;,(100 GeV) ~ 2.92 GeV versus mb^rrib) ~ 4.62 GeV yields (before cuts) a cross section of 
~ 0.75 fb versus ~ 2 fb, respectively, at tan /? = 50. The signal rates and resulting errors quoted 
in this section are those computed using = 4.62 GeV. Use of the running mass would reduce 
the event rates and increase our error estimates. In subsequent sections and figures, all results 
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FIG. 3: The tan P statistical error for C = 2000 fb"^ and tua = 100, 150, 200 GeV for 10% selection 
efficiency. For = 100 GeV, the signal and background rates are four times those given in Fig. |^ at 
the 10% efficiency point. Similar results are employed at = 150 and 200 GeV. 

and errors are computed in the MSSM context using the running 6-quark mass. Higher-order 
corrections of all kinds will be better known by the time the Linear Collider (LC) is constructed 
and data is taken and thus should not be a significant source of systematic uncertainty. An 
experimental challenge is associated with knowing the exact efficiency of the event selection 
procedure. At the working point of an efficiency esei = 10%, to achieve Atan/5/ tan/3 < 0.05 
requires Aegei/csei < 0.1, equivalent to Acgei < 1%- 

In addition to the hA Higgs boson background, two other Higgs boson processes could lead 
to a bbbb topology. First, the process e~^e~ —>■ HZ can give a bbbb final state. In fact, for large 
tan (3 the HZ cross section is maximal and similar in size to the hA cross section. Nonetheless, 
its contribution to the background is much smaller because the HZ — > bbbb branching is below 
10% compared to about 80% for hA bbbb. Since the hA process contributed only 2% 
of the total background, the contribution to the background from the HZ process can be 
neglected. The second Higgs boson process leading to a bbbb topology is that already discussed, 
e~^e~ — > bbh. The only distinction between this and the e~^e~ bbA process is a small difference 
in the angular distribution due to the different production matrix elements. Thus, the selection 
efficiency is almost identical. The production rate of the bbA process is proportional to tan^ /3 
while the bbh production rate is proportional to sin^ a/ cos^ /3. In the MSSM context, this 
latter factor is ~ tan^/3 for > 100 GeV and large tan/3 (assuming Msusy ~ 1 TeV). In 
the general 2HDM, since tan/3 ?a 1/ cos/3 at large tan/3, the expected rate depends mostly on 
sin a and the h mass. In this more general case where ^ but the MSSM expectation 
of a ~ — /3 ~ —tt/2 does not hold, the enhancement of the bbA signal by the bbh addition 
would only allow a determination of | sina| as a function of the presumed value of tan/3 (using 
the constraint that one must obtain the observed number of bbh + bbA events). Independent 
measurements of the HZ and hA production rates would then be needed to determine the value 
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of /3 — a and only then could a and (3 be measured separately. 

It is essential for the tan/3 determination that a very high integrated luminosity can be 
accumulated (we assume £ = 2000 fb"^ after several years of data-taking). Fig. ^ shows the 
expected statistical error on tan /? for = 100, 150 and 200 GeV, assuming that the only 
measured process is bbA with the help of our selection cuts. At the two higher values, in 
the MSSM context it is the H that would be decoupled and have mass rriH ~ ttlj^ and the h 
would be SM-like. Consequently, the bbH rate would be essentially identical to the bbA rate 
and, assuming that one could verify the MSSM Higgs context by independent means, would 
lead to still smaller tan (3 statistical errors than plotted, the exact decrease depending upon the 
signal-to-background ratio. For rriA = 150 and 200 GeV, the HA process (like the hA process 
at rriA = 100 GeV) would have to be computed in a specific model context or its relative weight 
fitted by studying bbbb production in greater detail in order to minimize any systematic error 
from this source. Results obtained in the case of the MSSM will be given in the following 
section. 

III. HA PRODUCTION: DECAY BRANCHING RATIOS AND TOTAL WIDTHS 

The branching ratios for H, A and decay to various allowed modes vary rapidly with tan /3 
in the MSSM when tan/? is in the low to moderate range, roughly below 20. Consequently, if 
these branching ratios can be measured accurately, tan /? can be determined with good precision 
in this range. Measurement of the branching fractions is most easily accomplished using HA and 
H^H^ pair production. In particular, the pair production processes are nearly independent 
of tan/5 so that the rate in a given channel provides a fairly direct probe of the branching 
ratio for that channel. That tan f] could be accurately determined using Higgs branching ratios 
measured in pair production was first demonstrated in |3^. Refs. [^, |3^ consider a number 
of models for which SUSY decays of the H, A and H^ are kinematically allowed. It was found 
that by measuring all available ratios of branching ratios it was possible to determine tan jS 
to better (often much better) than 10% for tan (3 values ranging from 2 up to as high as 25 
to 30 for rriA in the 200-400 GeV range, assuming y/s = 1 TeV and an effective luminosity 
(defined as the total luminosity times the selection efficiency of the cuts required to isolate 
the pair production process) of CcS = 80 fb^^ (equivalent, for example, to C = 2000 fb^^ 
for a selection efficiency of 4%). A more recent analysis using a few specific points in MSSM 
parameter space, focusing on the bbbb event rate and including a study at ^/s = 500 GeV, is 



given in l^g. This latter study uses a selection efficiency of 13% and negligible background for 
detection of e~^e~ — * hA bbbb (relevant for rriA < 100 GeV) or e~^e~ —>■ HA —>■ bbbb (relevant 
for tua > 150 GeV) and finds small errors for tan f3 at lower tan (3 values. Both []33| , Q and 



36| assume MSSM scenarios in which there are significant decays of the A and H to pairs of 



SUSY particles, in particular neutralinos and charginos. These decays remain non-negligible 
up to fairly high tan/? values. As a result, the bb branching fractions of the A and H increase 
more markedly as tan/? increases than if SUSY decays are absent. Indeed, in the absence of 
SUSY decays, the bbbb rate asymptotes quickly to a fixed value as tan /? increases. As we shall 
see, this means that much smaller errors for the tan (3 determination using the HA bbbb rate 
are achieved if SUSY decays are present. 

We now examine the errors on tan (3 that could be achieved using Higgs pair production, 
following procedures related to those of 0, ^ , but using updated luminosity expectations 



and somewhat more realistic experimental assumptions and analysis techniques. We restrict 
the analysis to the process e~^e~ — >■ HA bbbb, ignoring possible additional sensitivity through 
ratios relative to other final states. With both Higgs bosons reconstructed in their bb final state 
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as two back-to-back clusters of similar mass, backgrounds are expected to be negligible. All 
the results of this section are obtained using version 2.0 of HDECAY for computing the 
branching ratios and total widths of the Higgs bosons. 

To understand the sensitivity to the presence of SUSY decays of the heavy Higgs bosons, 
two different MSSM scenarios are considered: 

(I) mA = 200 GeV, m~ = 1 TeV, /i = Ma = 250 GeV, 

rrir = mr = mr = mr = mr= 1 TeV, 

= Ah = 0, At = ji/ tan/9 + \/Qmj (maximal mixing); 

(II) niA = 200 GeV, m~ = 350 GeV, /i = 272 GeV, Ma = 120 GeV, 
mr =mr = 356 GeV, mr = 273 GeV, mr = 400 GeV, 

= 0,Ab = -672 GeV, At = -369 GeV. 

In scenario (I), SUSY decays of the H and A are kinematically forbidden. Scenario (II) is taken 
from ||36[ in which SUSY decays (mainly to XiXi) ^ire allowed. We will assume that appropriate 



event selection criteria can be found such that for an event selection efficiency of 10% there will 
be negligible background. The resulting HA bbbb event rates (per 2000 fb^^ of integrated 
luminosity) are plotted for ^/s = 500 GeV in Fig. ^ as a function of tan f3. The difference in 
the dependence of the event rates on tan/3 is apparent. In more detail: in scenario (I) the 
bbbb event rates, after 10% selection efficiency, are 8, 77, 464, 1762, and 1859 at tan/5 = 1, 2, 
3, 10, and 40, respectively. The corresponding event rates in scenario (II) are 1, 5, 34, 1415 
and 1842. These differing tan j3 dependencies imply significant sensitivity of the tan j3 errors 
to the scenario choice, with worse errors for scenario (I). Finally, we note that for tan/? > 2 
the above event numbers are such that backgrounds are indeed negligible after 10% efficient 
selection cuts; for tan/? ~ 1, backgrounds might become an issue. 

To determine the la statistical errors of the tan (3 determination, we compute, for each choice 



of tan /3, the la upper and lower bounds on the expected event number as N{bbbb) ± y N{bbbb). 
These upper and lower bounds are also shown in Fig. ^. The upper (lower) event rate numbers 
are required to be > 10 to set an upper (lower) tan/3 bound, respectively. Since the event 
number increases monotonically with tan/3 for both MSSM scenarios, we can then use the 
given MSSM model scenario to determine the tan/? value for which the number of events is 
equal to the la upper (lower) bound. These values define the la upper (lower) bound on tan/3, 
respectively. The resulting fractional upper and lower limit errors Atan/3/ tan/? are plotted 
for MSSM scenarios (I) and (II) in Fig. |^. This procedure assumes that other measurements of 
SUSY particle production at the LHC and the LC will have fixed the MSSM scenario. 

Let us discuss in more detail the tan /3 errors from the HA — * bbbb rate in scenario (II) as 
compared to scenario (I). From Fig. ^ we see that in scenario (I) once tan/3 reaches 10 to 12 
the bbbb rate will not change much if tan (3 is increased further since the H ^ bb and A ^ bb 
branching ratios approach constant values. In contrast, if tan/? is decreased the bbbb rate 
declines significantly as other decay channels come into play. Thus, meaningful lower bounds 
on tan /3 are retained out to relatively substantial tan /3 values whereas upper bounds on tan /3 
disappear for tan /? > 10 — 12. In scenario (II), for reasons explained below, we have not plotted 
upper bounds on tan/3 for tan/? > 30. In fact, our numerical results indicate that the upper 
bound on A tan/?/ tan/? decreases again as tan/? increases beyond 30. We have traced this to 
the fact that HDECAY predicts that m^ decreases (at fixed m^ = 200 GeV) as tan j3 increases 
beyond 30. This results in an increase of the HA production cross section with increasing 
tan /?. This, in turn, implies that the bbbb rate increases (as shown in Fig. ^ and that we can 
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FIG. 4: Solid curves in the upper two figures give the rates for e~^e~ HA — > bbbb in scenarios 
(I) and (II). The solid curve in the lower figure is the resolved width F^^, Eq. (Q), for scenario (I). 
Dashed curves in all three figures correspond to the la upper and lower bounds on these quantities. 
We take niA = 200 GeV, ^/s = 500 GeV and C = 2000 fb~^ An efficiency of 10% is assumed 
for cuts, acceptance and tagging. The upper and lower la bounds for include an additional 
efficiency factor of 0.75 for keeping only events in the central mass peak and assume the estimated 
mass resolution of Pros = 5 GeV, including 10% systematic uncertainty. Results for ^ in SUSY 
scenario case (II) are very similar to those plotted for case (I). HDECAY [^] is used to compute the 
H and A widths and branching ratios. 
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FIG. 5: For the MSSM with niA = 200 GeV, and assuming £ = 2000 fb"^ at ^ = 500 GeV, we plot 
the la statistical upper and lower bounds, Atan/3/ tan/?, as a function of tan/3 based on: the rate 
for e+e" —>■ bbA + bbH bbbb; the rate for e^e~ HA — > bbbb; and the average resolved width 
defined in Eq. ^ for the H and A as determined in e^e~ — > HA bbbb events. For the rates, results 
for SUSY scenarios (I) and (II) differ significantly, as shown. For bbA + bbH — > bbbb and F^ ^ we show 
only the results for MSSM scenario (I). Results for scenario (II) are essentially identical. Upper and 
lower curves of a given type give the upper and lower la bounds, respectively, obtained using a given 
process as shown in the figure legend. We include running 6-quark mass effects and employ HDECAY 

H- 

obtain an upper bound on tan/5 despite the fact that the HA —>■ bbbb final state branching 
ratio approaches a constant value. However, since this predicted decrease of niH at high tan (3 
is somewhat peculiar to the precise parameters chosen for scenario (II), we do not regard this 
result as representative. For this reason, we have chosen not to show the scenario (II) upper 
limit curve beyond tan (3 = 30. Had we plotted the region above tan jS = 30, one would see a 
slowly declining upper limit on tan (3. 

The above results can be compared to the tan (3 determination based on the bbH+bbA bbbb 
rate using the procedures of Sec. II applied in the MSSM model context. For the computation 
of this rate, our calculation of the bbH and bbA cross sections includes the dominant radiative 
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corrections as incorporated via 6-quark mass running starting with mh{mb) = 4.62 GeV. The 
H and A branching ratios and widths are computed using HDECAY. Since there is httle 
sensitivity of this rate to the MSSM scenario (for the high tan f3 values for which this means 
of determining tan/3 is useful) we only present results for scenario (I); where plotted, errors for 
tan (3 from the bbH + bbA bbbb rate are essentially independent of the MSSM scenario choice. 
The errors on tan (3 resulting from the rate for bbH + bbA — > bbbb quickly become far smaller 
than those based on HA bbbb once tan /3 > 30. This is illustrated in Fig. ^, which compares 
the results for A tan /5/ tan j3 obtained using the e~^e~ —>■ HA — * bbbb rate to those based on the 
bbH + bbA — ^ bbbb rate. This comparison shows the natural complementarity between these two 
techniques for measuring tan/3. However, with these two techniques alone, there is always a 
range of intermediate-size tan j3 values for which a good determination of tan (3 is not possible. 

This "gap" can be partly filled, and the error on tan /5 at high tan /5 can be greatly reduced, 
by using the intrinsic total widths of the H and A to determine tan (3. However, it is only for 
tan (3 > 10 that the intrinsic widths can provide a tan (3 determination. This is because (a) the 
widths are only > 5 GeV (the detector resolution discussed below) for tan (3 > 10 and (b) the 
number of events in the bbbb final state becomes maximal once tan/5 > 10. 

We now discuss the experimental issues in determining the Higgs boson width. The expected 
precision of the SM Higgs boson width determination at the LHC and at a LC was studied 
in As described in a simple estimate (based on a detector energy flow resolution of 



res 



AE/E = 0.3/ vE for each of the two 6-jets) yields an expected detector resolution of Fr 
5 GeV for rriA ~ 200 GeV. However, an overall fit to the bb mass distribution similar to the one 
in the study of ||5^ would give a Higgs boson resonance peak width which is about 2a larger 



than that expected from the convolution of the 5 GeV resolution with the intrinsic Higgs width. 
This can be traced to the fact that the overall fit includes wings of the mass distribution that 
are present due to wrong pairings of the 6-jets. The mass distribution contains about 400 di-jet 
masses (2 entries per HA event), of which about 300 are in a central peak. If one fits only 
the central peak, the width is close to that expected based on simply convoluting the 5 GeV 
resolution with the intrinsic Higgs width. This indicates that about 25% of the time wrong jet- 
pairings are made and contribute to the wings of the mass distribution. Therefore, our estimates 
of the error on the determination of the Higgs width will be based on the assumption that only 
3/4 of the events {i.e. those in the central peak) retained after our basic event selection cuts 
(with assumed selection efficiency of 10%) can be used in the statistics computation. The m^,^ 
for each of the bb pairs identified with the H or A is binned in a single mass distribution. This is 
appropriate since the H and A are highly degenerate for the large tan (3 values being considered 
so that the resolution of 5 GeV is typically substantially larger than the mass splitting. Our 
effective observable is then the resolved average width defined by: 



pR _ 1 
^ H,A — 2 



res I 



(7) 



The resolved average width, F^^, for SUSY scenario (I) (including rriA = 200 GeV) is plotted 
in Fig. I as a function of tan (3. The results for scenario (II) are indistinguishable. 

In order to extract the implied tan (3 bounds, we must account for the fact that the de- 
tector resolution will not be precisely determined. There will be a certain level of systematic 
uncertainty which we have estimated at 10% of Fres, i-e. AF^^^ = 0.5 GeV. This systematic 
uncertainty considerably weakens our ability to determine tan/3 at the lower values of tan/3 
for which F^^. and F^^ are smaller than Fres- This systematic uncertainty should be carefully 
studied as part of any eventual experimental analysis. Given Fres, ^fres ^^"^ number of 
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selected bbbb events, N{bbbb), we compute the useful number of entries in the bbbb mass distri- 
bution for determining as A^entries = 2 x 0.75 x N{bbbb). The factor of 2 is because each 
bbbb event results in two entries, one for the H and one for the A, and the factor of 0.75 is that 
for retaining only the central peak of the distribution. The error for is then computed 
(following the procedure of |^) as 



H,A 



entries , 



1/2 



(8) 



The resulting upper and lower la bounds on F^_^ are plotted in Fig. |[ The upper and lower 
limits on tan f3 are then obtained as the values tan /3 ± A tan /3 for which the central prediction 
(the solid curve of Fig. |) agrees with the values F^^ ± AF^^. In computing AF^^ we have 
assumed a selection efficiency of 10% for computing N(bbbb). These errors are for C = 2000 fb~^ 
and a/s = 500 GeV. That an excellent determination of tan f3 will be possible at high tan (3 
is apparent. The resulting accuracy for tan/3 obtained from measuring the average (resolved) 
H/A width is shown in Fig. ^. We see that good accuracy is already achieved for tan/? as low 
as 25 with extraordinary accuracy predicted for very large tan/3. The sharp deterioration in 
the lower bound on tan /? for tan /? < 24 occurs because the width falls below F^es as tan j3 is 
taken below the input value and sensitivity to tan/? is lost. If there were no systematic error 
in Fres, this sharp fall off would occur instead at tan (3 < 14. To understand these effects in 
a bit more detail, we again give some numbers for scenario (II). At tan/3 = 50, 55 and 60, 
(F^(.,F;^^.) ~ 10.4, 12.5 and 14.9 GeV, respectively. After including the detector resolution, the 
effective average widths become 11.5, 13.4 and 15.7 GeV, respectively, whereas the total error 
in the measurement of the average width, including systematic error, is ~ 0.54 GeV. Therefore, 
tan/? can be determined to about ±1, or to better than ±2%. This high-tan/? situation can be 
contrasted with tan/? = 15 and 20, for which (F^^.,Fj^^) = 0.935 and 1.64 GeV, respectively, 
which become 5.09 and 5.26 GeV after including detector resolution. Meanwhile, the total 
error, including the statistical error and the systematic uncertainty for Fj-es, is about 0.57 GeV 
and no sensitivity to tan j3 is obtained. 

The accuracies from the width measurement are somewhat better than those achieved using 
the bbA + bbH bbbb rate measurement. However, both of these high-tan/? methods for 
determining tan /? are important because they are beautifully complementary in that they rely 
on very different experimental observables. Further, both methods are nicely complementary in 
their tan /? coverage to the tan (3 determination based on the HA bbbb rate, which comes in at 
lower tan/?. In fact, the width measurement can provide a decent tan/3 determination even in 
the previously identified "gap" region where neither the HA bbbb nor the bbA + bbH — > bbbb 
rates were able to provide such a determination. In particular, in the case of MSSM scenario 
(II), combining the HA bbbb rate and the width measurements implies that the worst tan/? 
lower bound is A tan /3/ tan/? ~ —0.25 at tan/3 ~ 28 and the worst tan/? upper bound is 
A tan/?/ tan /? ~ +0.30 at tan /? ~ 20. However, in the case of MSSM scenario (I) a good upper 
bound on tan /? is not possible if tan /? ~ 12 — 15, even after including the width measurement. 
Overall, there is a window, 10 < tan/? < 25 in scenario (I) or 20 < tan/? < 25 in scenario (II), 
for which an accurate determination of tan/? (A tan/?/ tan/? < 0.2) using just the bbbb final 
state processes will not be possible. This window expands rapidly as rriA increases (keeping s/s 
fixed). Indeed, as increases above 250 GeV, HA pair production becomes kinematically 
forbidden at -y/s = 500 GeV and, in addition, detection of the bbH + bbA processes at the LC 
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(or the LHC) requires [0] increasingly large values of tan/5. This difficulty persists even for 
•y/s ~ 1 TeV and above; if ttia > \fsl'^i the H and A cannot be pair-produced and yet the rate 
for U)H + 66^4 production is undetectably small for moderate tan /? values. 

In the above study, we have not made use of other decay channels of the H and A, such 



as if WW, ZZ, H ^ hh, A ^ Zh and H,A ^ SUSY. The theoretical studies of 34 



indicate that their inclusion could improve the precision with which tan/3 is measured at low 
to moderate tan P values. A determination of ^ is also possible using the bbA + bbH — > bbbb 
events. To estimate how well tan/3 can be determined in this way, let us assume that 50% 
of the events selected in the analysis of Sec. II can be used for a fit of the average width and 
that (as in the HA bbbb study) a resolution of 5 GeV can be achieved, based on a detector 
energy fiow resolution of AE/E = O.S/Ve. If we again assume that 10% systematic error for 
the width measurement can be achieved, the resulting tan/? errors would be similar to those 
obtained from the bbA + bbH — » bbbb event rate for tan/3 > 30 (see Fig. |^), i.e. not as small 
as the tan/3 errors obtained from the measurement of F^^ in HA bbbb events. A complete 
analysis that takes into account the significant background and the broad energy spectrum of 
the radiated H and A is needed to reliably assess the tan/3 errors that would be obtained by 
measuring F^ ^ in bbA + bbH — > bbbb events. If tan /3 is large enough, bbA + bbH bbbb would 
be observable in the MSSM even when > \^/2, or, for a more general model, whenever 
HA pair production is kinematically forbidden but bbA — »• bbbb and/or bbH bbbb production 
is allowed. Then, the event rate for, and width measurements from bbA and/or bbH production 
would allow a determination of tan/?. We have not attempted a quantitative study of this 
situation. 

Let us briefiy return to the interpretation of these measurements in terms of tan /?. As stated 
in the introduction, we are using tan /? as a tree-level mnemonic to characterize the bb Yukawa 
coupling of the Higgs bosons. For the soft-SUSY-breaking parameters for MSSM scenarios (I) 
and (II), the one- loop corrections to the bb couplings of the H and A and the stop/sbottom 



mixing present in the one-loop corrections to the Higgs mass matrix |Q are small. More 
generally, however, substantial ambiguity can arise, especially if the sign and magnitude of fi 
is not fixed. Assuming that these parameters are known, the errors for the Yukawa coupling 
obtained from these measurements can be related to any given definition of tan /? and, except 
in very unusual cases, the resulting error on tan/3 would be fairly insensitive to the precise 
scenario. For example, one possible definition of tan/3 would be that the fi~^fi~A coupling 
should be precisely given by —(rrife/f) tan /3, see Eq (H). This is a convenient definition since 
the ^~^fi~A coupling will have very modest higher-order corrections relative to the tree- level 
and any such corrections can then be sensibly absorbed using the above definition of tan/3. 
Given this definition of tan /3, the Hbb and Abb couplings can be computed to any desired order 
once the necessary MSSM parameters are known. In this way, all the probes of heavy Higgs 
Yukawa couplings discussed in this paper can be related to this common definition of tan /?. 



IV. H+H- PRODUCTION: DECAY BRANCHING RATIOS AND TOTAL WIDTH 

In this section, we extend our study to include charged Higgs boson production processes. 
Existing analyses of e~^e^ —>■ H^H^ production indicate that the absolute event rates and 
ratios of branching ratios in various H^H~ final state channels will allow a relatively accurate 
determination of tan/? at low tan/3 ||33|, 0. The process e^e~ —>■ H'^tb can also be sensitive 



to tan/3 [35|. Here, we focus on an experimentally based analysis of the determination of tan/? 
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using the H^H~ tbtb event rate. As anticipated on the basis of the earher work referenced 
above, we find that good accuracy can be achieved at low tan/5. We also demonstrate that 
the total width of the measured in the tb decay channel using H^H~ — > tbtb production 
will allow a fairly precise determination of tan (3 at high tan (3. Since these two techniques for 
determining tan (3 are statistically independent of one another and of the tan j3 measurements 
that employ neutral Higgs production, they will increase the overall accuracy with which tan j3 
can be measured at both low and high tan (3. 

The reaction e^e'^ — > H^H~ tbtb can be observed at a LC [5^, and recent high- 



luminosity simulations |Q show that the cross section times branching ratio can be measured 
precisely. As soon as the charged Higgs boson decay into tb is allowed this decay mode is dom- 
inant. Nonetheless, BR(if^ — >■ tb) varies significantly with tan/3, especially for small values of 
tan 13 where the tb mode competes with the rz/ mode. The tb branching ratio and width 

are sensitive to tan f3 in the form 

r(if± ^ t6) oc m^cot^/J + m^tan^/?. (9) 

As in Sec. Ill, we will use HDECAY (which incorporates running of the 6-quark mass) to 
evaluate the charged Higgs boson branching ratios and decay width. It is useful to note that 
the above form results in a minimum in the tb partial width and branching ratio in the vicinity 
of tan/3 ~ 6 — 8. The depth of the minimum in the branching ratio depends upon the extent 
to which the tb channel is competing against other modes. In contrast, the cross section for 



e^e^ H^H^ production is independent of tan [3 at tree-level. (The one-loop corrections I^T 
result in a 10% variation of the cross section with tan/3 which must be taken into account 
when the data are taken; we do not include them in our study.) The net result is that the rate 
for e"'"e~ H^H~ — > tbtb has significant dependence upon tan/3, coming mainly from the 
variation in the branching ratio. 

Our procedures for estimating errors for the tbtb rate and for the total width are similar 
to those given earlier for HA production rates and width in the hbhb channel. We base our 



efficiency for the tbtb final state on the study of [Q. For mH± = 300 GeV at y/s = 800 GeV, 
this study finds that the tbtb final state can be isolated with an efficiency of 2.2%. The reason 
for the much lower efficiency as compared to 10% efficiency for the bbbb final state of HA 
production is the difficulty of assigning the non-6 jets from t-decays to the correct top mass 
cluster. For = 500 GeV and mH± = 200 GeV, we have adopted the same 2.2% efficiency, 
for which we assume little or no background after cuts. For the total width determination, 
we assume that we keep only 75% of the events after cuts {i.e. a fraction 0.75 x 0.022 of the 
raw event number), corresponding to throwing away wings to the mass peaks, and each tbtb 
event is counted twice since we can look at both the H~^ and the H~ decay. We define a 
resolved width which incorporates the intrinsic resolution for the width determination, taken 
to be Frp, = 5 GeV: 



r'^^ = \jnt?+[rres?. (10) 

Estimated errors based on the width measurement will assume a 10% systematic error in our 
knowledge of Fres, i-e. AF^^| = 0.5 GeV as for the H,A case. We employ Eq. (||, with the 
replacement F^^ — > F^±, to compute AF^±. In this case, A'entries = 2 x 0.75 x N{tbtb), where 
N{tbtb) is computed using the above-noted selection efficiency of 0.022. Figure ^ shows the 
resulting tbtb final state rate, N{tbtb), for MSSM scenarios (I) and (II) and the resolved width 
(F^±) for scenario (I). Also shown are the corresponding la upper and lower bounds on the 
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FIG. 6: The solid curves give the rates for e+e^ — > H^H^ tbtb for MSSM scenarios (I) and (II) 
in the upper two figures and r|±, Eq. (^), for scenario (I) in the lower figure. The dashed curves 
are the corresponding la upper and lower bounds. We take infj± = 200 GeV, ^/s = 500 GeV and 
C = 2000 fb"^ An efficiency of 2 .2% is assumed for cuts, acceptance and tagging. The upper and 
lower la bounds for include an additional efficiency factor of 0.75 (which corresponds to keeping 
only events in the central mass peak) and assume the estimated mass resolution of Fres = 5 GeV, 
including 10% systematic uncertainty. The F^± results obtained in scenario (II) are very similar to 
those plotted for scenario (I). 
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rate and resolved width. These are then used in exactly the same manner as described in the 
HA case to determine the upper and lower bounds on tan p. 

For the rate, we observe from Fig. |^ that upper bounds on tan f3 will be poor once tan /? > 10 
because of the very slow variation of the tbtb final rate in this region. At high tan/?, for SUSY 
scenario (1) lower bounds will be determined by the part of the tbtb rate curve that rises 
rapidly when tan/? falls below 10. For SUSY scenario (II), the beginning of the dip will fix the 
lower bounds on tan (3 when tan/3 is large, assuming that we know ahead of time from other 
experimental data that tan/5 is larger than 15. As regards the width, the main point to note 
is that r^± rises only slowly with increasing tan/5. As a result, the 5 GeV resolution and 10% 
systematic error for this resolution are significant compared to the < 10 GeV width that 
applies throughout the tan/3 range studied. Note also, that for moderate tan/5 values, there 
will be no lower bound on tan /5 as a result of the fact that r^± never falls below Fres = 5 GeV, 
while the la errors are substantially lower than this. We will also assume that if tan f3 is large, 
then we will know from other experimental information (such as the HA final state) that tan (3 
is not small and that the small rise in the width for tan/5 ~ 1 is not relevant. 

The resulting tan/5 upper and lower bounds appear in Fig. |^. Comparing to Fig. ^, we 
observe that for SUSY scenario (I) the tbtb rate measurement gives a tan/5 determination that 
is quite competitive with that from HA production in the bbbb final state. For SUSY scenario 
(II), the tbtb rate gives an even better tan/5 determination than does the bbbb rate. On the 
other hand, the width measurement from the tbtb final state of H^H~ production is much 
poorer than that from the hbbb final state of HA production, as was to be expected from the 
discussion given earlier. 

The rate for e"'"e~ tbH~ + tbH^ — > tthb is also very sensitive to tan (3 and might be 
a valuable addition to the e"'"e~ H^H~ — > tbtb and e'^e~ bbA + bbH — >• bbbb rate 
determinations of tan (5. The theoretical study of I^S] finds, for example, that if = 
200 GeV and tan/5 = 50 (tan/5 = 20), then the la errors (including systematic uncertainties) 
on tan/5 are Atan/5/ tan/5 = 0.06 (Atan/5/ tan/5 = 0.2), respectively, for C = 2000 fb""*^ and 
^ = 500 GeV. 

By combining in quadrature the tan/5 errors for the various individual measurements, as 
given in Figs. ^ and 0, we obtain the net errors on tan j3 shown in Fig. ^. 

V. COMPARISON TO LHC DETERMINATIONS OF tan (3 

In this section, we will compare the LC results summarized in Figs. ^ ^ and |^ to the tan /5 
accuracies that can be achieved at the LHC based on H, A, H^ production and decay processes. 
First note that there is a wedge-shaped window of moderate tan (3 and > 200 GeV for which 
the A, H and if^ are all unobservable (see, for example, Refs. ||21], ^ ^]). In this wedge, the 
only Higgs boson that is detectable at the LHC is the light SM-like Higgs boson, h. Precision 
measurements of the properties of the h typically only provide weak sensitivity to tan/5, and 
will not be considered here. The lower tan/5 bound of this moderate-tan /5 wedge is defined 
by the LEP-2 limits [Q, which are at tan/5 ~ 3 for itla ~ 200 GeV, falling to tan/5 ~ 2.5 for 
rriA ^ 250 GeV, assuming the maximal mixing scenario [see SUSY scenario (I) defined earlier] . 
The upper tan /? limit of the wedge is at tan /5 ~ 7 for mA ~ 200 GeV rising to tan /5 ~ 15 at 
TJiA ~ 500 GeV. For either smaller or larger tan /? values, the heavy MSSM Higgs bosons can 
be detected and their production rates and properties will provide sensitivity to tan /?. 

We will now summarize the results currently available regarding the determination of tan /5 
at the LHC using Higgs measurements (outside the wedge region) assuming a luminosity of 
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FIG. 7: For the MSSM with mjy± ~ = 200 GeV, and assuming C = 2000 fb"^ at = 500 GeV, 
we plot the Icr statistical upper and lower bounds, A tan (3/ tan /3, as a function of tan (3 based on: the 
rate for e~^e~ H^H^ — > tbtb; and the resolved width F^j- defined in Eq. ( |To| ) as determined in 
e^e~ — > H^H^ tbtb events. For the rates, results for SUSY scenarios (I) and (II) differ significantly, 
as shown. For we show only the results for MSSM scenario (I). Results for scenario (II) are 
essentially identical. Upper and lower curves of a given type give the upper and lower la bounds, 
respectively, obtained using a given process as shown in the figure legend. We include running 5-quark 
mass effects and employ HDECAY 



C = 300 fb . The methods employed are those proposed in ^0 
studied at the LHC are the following. 



The reactions that have been 



H ^ ZZ ^Ai^. 

The best accuracy that can be achieved at low tan (3 is obtained from the H ZZ A£ 
rate. One finds A tan/?/ tan/5 = ±0.1 at tan/3 = 1 rising to > ±0.3 by tan/? = 1.5 for 
the sample choice of mn = 300 GeV. For rriH < '^^rriz, tan/5 cannot be measured via this 
process. In the MSSM maximal mixing scenario, such low values of tan (3 are unlikely in 
light of LEP-2 results. 



2. gg^H + gg^ A 



T ' T 



, fi^fi and gg bbH ± bbA —>■ bbr^r , bbfi^fi |22 
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FIG. 8: For the MSSM with mj^± ~ = 200 GeV, and assuming C = 2000 fb"' at = 500 GeV, 
we plot the la statistical upper and lower bounds, Atan/3/ tan/3, as a function of tan/3 based on 
combining (in quadrature) the results shown in Figs. |5| and |^. Results are shown for the SUSY 
scenarios (I) and (II) described in the text. 



At high tan/? and taking rriA = 150 GeV, Fig. 19-86 of shows that the gg ^ H ^ 
T^T~ , gg —>■ A T^T~ , and gg — ^ bbA + bbH —>■ bbT^r" rates can, in combination, be 
used to determine tan/5 with an accuracy of ±0.15 at tan/3 = 5, improving to ±0.06 
at tan/5 = 40. The corresponding rates with H,A —>■ fi'^fi~ yield a somewhat better 
determination at higher tan/3: ±0.12 at tan/3 = 10 and ±0.05 at tan/3 = 40. 

Interpolating, using Figs. 19-86 and 19-87 from [^, we estimate that at ttia ~ 200 GeV 
(our choice for this study) the error on tan f3 based on these rates would be smaller than 
±0.1 for tan/3 > 13, asymptoting to ±0.05 at large tan/3. 

It is important to note that the tan (3 sensitivity for tan /3 < 20 — 30 is largely due to the 
loop-induced gg —>■ A and gg ^ H production processes. Thus, interpreting these fully 
inclusive rates in terms of tan /3 when tan /3 < 20 — 30 requires significant knowledge of 
the particles, including SUSY particles, that go into the loops responsible for the gg ^ H 
and gg A couplings. 
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The importance of including the gg ^ H, gg ^ A as well as the gg — * bbH + gg ^ bbA 
processes in order to obtain observable signals for tan j3 values as low as 10 in the fi~^fi~ 
channels is apparent from [Q. For C = 300 fb~^ and tua = 200 GeV, they find that the 
bbfj,~^fj,~ final states can only be isolated for tan/3 > 30 whereas the inclusive //"""/i" final 
state from all production processes becomes detectable once tan/3 > 10. 



3. gg tbH- + tbH+ with r±z/ [^|. 

The tbH^ tbru rate gives a fractional tan/5 uncertainty, Atan/5/ tan/3, ranging from 
±0.074 at tan/3 = 20 to ±0.054 at tan/? = 50. This signal is somewhat cleaner to 
interpret in terms of a tan/3 measurement than the inclusive signals for the H and A 
summarized above, since there are no uncertainties related to SUSY loop contributions. 

Sensitivity to tan (3 deriving from direct measurements of the decay widths has not been studied 
by the LHC experiments. One can expect excellent tan/3 accuracy at the higher tan/? values 
for which the gg — > bbfi^fi~ signal for the H and A is detectable. 

Let us now compare these LHC results to the LC errors for tan /3, assuming = 200 GeV. 
First, consider tan /3 < 10. As summarized above, the LHC error on tan j3 is ±0.12 at tan /3 ~ 10 
and at tan/3 ~ 1, and the error becomes very large for 1.5 < tan/3 < 5. Meanwhile, the LC 
error from Fig. |^ ranges from roughly ±0.03 to ±0.05 for 2 < tan/3 < 5 rising to about 
±0.1 at tan/3 ~ 10 [in the less favorable SUSY scenario (I)]. Therefore, for tan/3 < 10 the 
LC provides the best determination of tan/3 using Higgs observables related to their Yukawa 
couplings. (In the MSSM context, other non-Higgs LHC measurements would allow a good 
tan f5 determination at low to moderate tan (3 based on other kinds of couplings.) In the middle 
range of tan/? (roughly 13 < tan/3 < 30 at rriA ~ 200 GeV), the heavy Higgs determination 
of tan/? at the LHC might be superior to that obtained at the LC. This depends upon the 
SUSY scenario: if the heavy Higgs bosons can decay to SUSY particles, the LC will give 
tan /? errors that are quite similar to those obtained at the LHC; if the heavy Higgs bosons do 
not have substantial SUSY decays, then the expected LC tan/3 errors are substantially larger 
than those predicted for the LHC. At large tan/?, the LC measurement of the heavy Higgs 
couplings and the resulting tan/? determination at the LC is numerically only slightly more 
accurate than that obtained at the LHC. For example, both are of order ±0.05 at tan/? = 40. 



TABLE I: A comparison of fractional errors, A tan/?/ tan/?, achievable for £ = 2000 fb"^ at the LC 
with those expected at the LHC for £ = 300 fb~^, assuming niA = 200 GeV in the MSSM. LC results 
are given for both SUSY scenarios (I) and (II), where Higgs boson decays to SUSY particles are 
disallowed, respectively allowed. LHC results are estimated by roughly combining the determinations 
of tan/3 based on H, A production from |22| with those using production from |24], both of which 
assume the standard MSSM maximal mixing scenario. All entries are approximate. 



tan /? range 


LHC 


LC (case I) 


LC (case II) 


1 


0.12 


0.15 


0.1 


1.5-5 


very large 


0.03-0.05 


0.03-0.05 


10 


0.12 


0.1 


0.05 


13-30 


0.05 


0.6-0.1 


0.05-0.1 


40-60 


0.05-0.03 


0.05-0.025 


0.05-0.025 
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These comparisons are summarized in Table |. It is possible that the net LHC tan/3 error 
would be somewhat smaller than the LC error for tan /3 > 40 if both ATLAS and CMS can 
each accumulate £ = 300 fb~^ of luminosity; combining the two data sets would presumably 
roughly double the statistics and decrease errors by a factor of order 1/ \/2. In any case, a very 
small error on tan (3 will be achievable for all tan (3 by combining the results from the LC with 
those from the LHC. |6| 

VI. CONCLUSIONS 

A high-luminosity linear collider will provide a precise determination of the value of tan/3 
throughout much of the large range of possible interest, 1 < tan/5 < 60. In this paper, we 
have studied the sensitivity to tan (3 that will result from measurements of heavy Higgs boson 
production processes, branching fractions and decay widths. These are all directly determined 
by the ratio of vacuum expectation values that defines tan /5, and each can be very accurately 
measured at an LC over a substantial range of relevant tan/5 values. In particular, there are 
several Higgs boson observables which are likely to provide the most precise measurement of 
tan/3 when tan/? is very large. In the context of the MSSM, there is a particularly large 
variety of complementary methods that will allow an accurate determination of tan/3 when 
^ so that e+e^ HA pair production is kinematically allowed. Using the sample 

case of rriA = 200 GeV (in the MSSM context) and a LC with y/s = 500 GeV, we have 
demonstrated the complementarity of employing: 

a) the bbA, bbH bbbb rate; 

b) the HA bbbb rate; 

c) a measurement of the average H, A total width in HA production; 

d) the H^H — > tbtb rate; and 

e) the total H^ width measured in H^H^ —>■ tbtb production. 

By combining the tan /3 errors from all these processes in quadrature, we obtain the net errors 
on tan j3 shown in Fig. || by the lines [solid for SUSY scenario (I) and dashed for SUSY scenario 
(II)], assuming a multi-year integrated luminosity of C = 2000 fb~\ We see that, independent 
of the scenario, the Higgs sector will provide an excellent determination of tan j3 at small and 
large tan /? values, leading to an error on tan /? of 10% or better. If SUSY decays of the H, A, H^ 
are significant [SUSY scenario (H)], the tan/? error will be smaller than 13% even in the more 
difficult moderate tan/? range. However, if SUSY decays are not significant [SUSY scenario (I)] 
there is a limited range of moderate tan /3 for which the error on tan (3 would be large, reaching 
about 50%. 

In the preceding section, we considered how these tan/3 errors from the LC compared to 
tan 13 errors determined at the LHC based only on measurements involving H, A, H^ production 
and decay. The broad conclusions were: (i) for low tan/? (< 10) the errors on tan/? from LHC 
Higgs measurements would be much larger than those attainable at the LC; (ii) for high tan j3 
(> 30) the LHC and LC tan/? errors were both small and quite comparable in magnitude; 
and (iii) in the moderate-tan /3 range (13 < tan/3 < 30) the LHC errors on tan/3 would very 
possibly be smaller than the LC errors. However, we also noted that in this latter region 
some care in interpretation of the LHC results would be necessary due to the need to include 
loop-induced gg H and gg ^ H production processes in order to obtain good sensitivity to 
tan/?; these might be infiuenced by loops of SUSY particles and, possibly, other undiscovered 
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new physics. This LHC versus LC comparison should also be viewed as highly preliminary 
since the LHC collaborations have not yet studied all the relevant observables. In particular, 
they have not looked at the tan (3 determination using the directly measured widths of the H 
and A. Regardless of the relative magnitude of the LHC versus LC tan (3 errors, the clean LC 
environment will provide an important and independent measurement that will complement any 
LHC determination of tan j3. Different uncertainties will be associated with the determination 
of tau/S at a hadron and an e^e~ collider because of the different backgrounds. Further, the 
LHC and LC measurements of tan j3 will be highly complementary in that the systematic errors 
involved will be very different. 

Combining all the different LC measurements as above does not fully account for the fact 
that the "effective" tan (3 value being measured in each process is only the same at tree-level. 
The tan /5 values measured via the if — > 66 Yukawa coupling, the A ^ bb Yukawa coupling 
and the — > tb Yukawa coupling could all be influenced differently by the MSSM one-loop 
corrections. For some choices of MSSM parameters, the impact of MSSM radiative corrections 
on interpreting these measurements can be substantial However, if the masses of the 

SUSY particles are known, so that the important MSSM parameters entering these radiative 
corrections (other than tan/3) are fairly well determined, then a uniform convention for the 
definition of tan/? can be adopted and, in general, an excellent determination of tan/? (with 
accuracy similar to that obtained via our tree-level procedures) will be possible using the linear 
collider observables considered here. Even for special SUSY parameter choices such that one 
of the Yukawa couplings happens to be significantly suppressed, the observables a)-e) would 
provide an excellent opportunity for pinning down all the Yukawa couplings and checking the 
consistency of the MSSM model. 

Finally, it is important to note that the techniques considered here can also be employed in 
the case of other Higgs sector models. For example, in the general (non-SUSY) 2HDM, if the 
only non-SM-like Higgs boson with mass below ^/s is the A , then a good determination of 
tan (3 will be possible at high tan (3 from the bbA bbbb production rate. Similarly, in models 
with more than two Higgs doublet and/or triplet representations, the Yukawa couplings of the 
Higgs bosons, and, therefore, the analogues of the 2HDM parameter tan/3, will probably be 
accurately determined through Higgs production observables in e~^e~ collisions. 
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